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Pion Physics at Low Energy and High Accuracy 

H. Leutwyler 

Institute for Theoretical Physics, University of Bern, Sidlerstr 5, CH-301 2 Switzerland 

Abstract. The role of the quark condensate for the low energy structure of QCD is discussed in 
some detail. In particular, the dependence of Mjt on m,, and nid and the low energy theorems for the 
KK scattering amplitude are reviewed. The new data on K^^ decay beautifully confirm the standard 
picture, according to which the quark condensate is the leading order parameter of the spontaneously 
broken chiral symmetry. 
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>• [ At energies that are small compared to {MwtMz,Mh} — O(lOOGeV), the weak 

^ ' interaction freezes out, because these energies do not suffice to excite the relevant 

(T^ ■ degrees of freedom. As a consequence, the gauge group of the Standard Model, 

C^: ! SU(3)xSU(2)xU(l), breaks down to the subgroup SU(3)xU(l) - only the photons, the 

gluons, the quarks and the charged leptons are active at low energies. Since the neutrini 

neither carry colour nor charge, they decouple. 

7h '. The effective Lagrangian relevant at low energies is the one of QCD + QED. The 

P^l strength of the interaction is characterized by the two coupling constants g and e. In 

p ■ contrast to the Standard Model, the SU(3)xU(l) Lagrangian does contain mass terms: 

(D . the quark and lepton mass matrices m^, m^. The field basis may be chosen such that m^ 

and m£ are diagonal and positive. 

The two gauge fields behave in a qualitatively different manner: while the photons do 
S^ . not carry electric charge, the gluons do carry colour. This difference is responsible for the 

fact that the strong interaction becomes strong at low energies, while the electromagnetic 
interaction becomes weak there, in fact remarkably weak: the photons and leptons 
essentially decouple from the quarks and gluons. For the QCD part, on the other hand, 
perturbation theory is useful only at high energies. In the low energy domain, the strong 
interaction is so strong that it confines the quarks and gluons. For the same reason, a term 
in the Lagrangian of the form ~ G^y G^^ (where G^y is the gluon field strength) cannot 
a priori be dismissed, despite the fact that it represents a total derivative: it generates an 
electric dipole moment in the neutron, for instance. Conversely, the experimental fact 
that the dipole moment is smaller than lO^^ecm implies that (in the basis where the 
quark mass matrix is diagonal, real and positive) the vacuum angle must be tiny, so 
that the Lagrangian is invariant under the discrete symmetries P, C and T, to a very high 
degree of accuracy. 

The resulting effective low energy theory is mathematically even more satisfactory 
than the Standard Model as such - it does not involve scalar degrees of freedom and 
has fewer free parameters. Remarkably, this simple theory must describe the structure 



of cold matter to a very high degree of precision, once the parameters in the Lagrangian 
are known. It in particular explains the size of the atoms in terms of the scale 



471 



e^nie 



which only contains the two parameters e and m^ - these are indeed known to an 
incredible precision. Unfortunately, our ability to solve the QCD part of the theory is 
rather limited - in particular, we are still far from being able to demonstrate on the basis 
of the QCD Lagrangian that the strong interaction actually confines colour. Likewise, our 
knowledge of the magnitude of the light quark masses leaves to be desired - we need to 
know these more accurately in order to test ideas that might lead to an understanding of 
the mass pattern, such as the relations with the lepton masses that emerge from attempts 
at unifying the electroweak and strong forces. 



2. SYMMETRIES OF MASSLESS QCD 

In the following, I focus on the QCD part and switch the electromagnetic interaction 
off. It so happens that the interactions of u,d,s with the Higgs fields are weak, so 
that the masses mu.nid.nis are small. Let me first set these parameters equal to zero 
and, moreover, send the masses of the heavy quarks, mc,mi,,mt to infinity. In this 
limit, the theory becomes a theoreticians paradise: the Lagrangian contains a single 
parameter, g. In fact, since the value of g depends on the running scale used, the 
theory does not contain any dimensionless parameter that would need to be adjusted to 
observation. In principle, this theory fully specifies all dimensionless observables as pure 
numbers, while dimensionful quantities like masses or cross sections can unambiguously 
be predicted in terms of the scale Aqcd or the mass of the proton. The resulting theory - 
QCD with three massless flavours - is among the most beautiful quantum field theories 
we have. I find it breathtaking that, at low energies, nature reduces to this beauty, as soon 
as the dressing with the electromagnetic interaction is removed and the Higgs condensate 
is replaced by one that does not hinder the light quarks, but is impenetrable for W and Z 
waves as well as for heavy quarks. 

The Lagrangian of the massless theory, i^g^D' has a high degree of symmetry that orig- 
inates in the fact that the interaction among the quarks and gluons is flavour-independent 
and conserves helicity: ig^^ is invariant under independent flavour rotations of the three 
right- and left-handed quark fields. These form the group G = SU(3)^ x SU(3)^. The 
corresponding 16 currents V^^qy^^Xiq and A^ =qy^y5j^i(l are conserved, so that their 
charges commute with the Hamiltonian: 

Vafa and Witten [1] have shown that the state of lowest energy is necessarily invariant 
under the vector charges: Q^ \0) = 0. For the axial charges, however, there are the two 
possibilities characterized in table 1. 

The observed spectrum does not contain parity doublets. In the case of the lightest 
meson, the pion, for instance, the lowest state with the same spin and flavour quantum 



TABLE 1. Alternative realizations of the symmetry group G = SU(3)^ x 811(3)^ 



e?|o)=o 


<2f|o)^o 


Wigner-Weyl realization of G 


Nambu-Goldstone realization of G 


ground state is symmetric 


ground state is asymmetric 


m,q,\o)^o 


{0\q,q,\0)^0 


ordinary symmetry 


spontaneously broken symmetry 


spectrum contains parity partners 


spectrum contains Goldstone bosons 


degenerate multiplets of G 


degenerate multiplets of SU(3)i, C G 



numbers, but opposite parity is the ao(980) . So, experiment rules out the first possibility: 
for dynamical reasons that yet remain to be understood, the state of lowest energy is an 
asymmetric state. Since the axial charges commute with the Hamiltonian, there must be 
eigenstates with the same energy as the ground state: 



^ecz,0|O)=ef//°cJO)=O. 



The spectrum must contain 8 states Qj |0), . . . , 2^ |0) with E = P = 0, describing mass- 
less particles, the Goldstone bosons of the spontaneously broken symmetry. Moreover, 
these must carry spin 0, negative parity and form an octet of SU(3). 



3. QUARK MASSES AS SYMMETRY BREAKING PARAMETERS 



Indeed, the 8 lightest hadrons, 'K^,'nP,n^,K^,K'^,K'^,K^ ,r\, do have these quantum 
numbers, but massless they are not. This has to do with the deplorable fact that we are 
not living in paradise: the masses mu.nid.nis are different from zero and thus allow the 
left-handed quarks to communicate with the right-handed ones. The full Hamilitonian is 
of the form 



H. 



QCD 



^QCD ' ^QCD 1 



Tjl 

^QCD 






The quark masses may be viewed as symmetry breaking parameters: the QCD- 
Hamiltonian is only approximately symmetric under independent rotations of the right- 
and left-handed quark fields, to the extent that these parameters are small. Chiral 
symmetry is thus broken in two ways: 



spontaneously 
explicitly 



{0\qRq,\0)^0 
mu^nicumsT^O 



The consequences of the fact that the explicit symmetry breaking is small may be worked 
out by means of an effective field theory, "chiral perturbation theory" [2, 3, 4]. In this 
context, the heavy quarks do not play an important role - as the corresponding fields 
are singlets under SU(3)rxSU(3)£, we may include their contributions in the symmetric 
part of the Hamiltonian, irrespective of the size of their mass. 



Since the masses of the two lightest quarks are particularly small, the Hamiltonian of 
QCD is almost exactly invariant under the subgroup SU(2)«xSU(2)£. The ground state 
spontaneously breaks that symmetry to the subgroup SU(2)v, - the good old isospin 
symmetry discovered in the thirties of the last century [5]. The pions represent the 
corresponding Goldstone bosons [6], while the kaons and the r\ remain massive if the 
limit niu^mj ^ is taken at fixed m^. In the following, I consider this framework and, 
moreover, ignore isospin breaking, setting niu = nid = m. 

If SU(2)« X SU(2)£ was an exact symmetry, the pions would be strictly massless. Ac- 
cording to Gell-Mann, Oakes and Renner [7], the square of the pion mass is proportional 
to the product of the quark masses and the quark condensate: 

Ml--^x{m^ + md)x\{0\Uu\0)\. (1) 

The factor of proportionality is given by the pion decay constant Fjt. The term m^ + m^ 
measures the explicit breaking of chiral symmetry, while the quark condensate, 

(0|mm|0) = {0\urUl\0)+c.c. = (0| JJ|0) , 

is a measure of the spontaneous symmetry breaking: it may be viewed as an order 
parameter and plays a role analogous to the spontaneous magnetization of a magnet. 



4. ROLE OF THE QUARK CONDENSATE 

The approximate validity of the relation (1) was put to question by Stem and collabora- 
tors [8], who pointed out that there is no experimental evidence for the quark condensate 
to be different from zero. Indeed, the dynamics of the ground state of QCD is not under- 
stood - it could resemble the one of an antiferromagnet, where, for dynamical reasons, 
the most natural candidate for an order parameter, the magnetization, happens to vanish. 
There are a number of theoretical reasons indicating that this scenario is unlikely: 

(i) The fact that the pseudoscalar meson octet satisfies the Gell-Mann-Okubo formula 
remarkably well would then be accidental. 

(ii) The value obtained for the quark condensate on the basis of QCD sum rules, in 
particular for the baryonic correlation functions [9], confirms the standard picture. 

(iii) The lattice values [10] for the ratio nis/m agree very well with the result of the 
standard chiral perturbation theory analysis [11], corroborating this picture further. 

Quite irrespective, however, of whether or not the scenario advocated by Stern et al. is 
theoretically appealing, the issue can be subject to experimental test. In fact, significant 
progress has recently been achieved in this direction [12, 13]. The remainder of the talk 
concerns this matter. 

The Gell-Mann-Oakes-Renner formula is not exact. The expansion of M^ in powers 
of mujinj contains an infinite sequence of contributions. The expansion starts with a 
term linear in the quark masses: 

Ml = M^-:^:^^M^ + 0{M^), M^ = {m, + md)B. (2) 



The coefficient B of the linear term is given by the value of | (0| uu \0)\/F^ in the limit 
i^u: f^d -^ 0, and F is the value of F^j in that limit. The Gell-Mann-Oakes-Renner formula 
is obtained by dropping the higher order contributions. These are dominated by the term 
of order M^, which involves one of the coupling constants occurring in the effective 
Lagrangian at order p^. More precisely, the formula involves the value of the running 
coupling constant £3 at scale n = M, which logarithmically depends on M. Expressed in 
terms of the corresponding intrinsic scale A3, we have 

h = ln§. (3) 

The symmetry does not determine the numerical value of this scale. The crude estimates 
underlying the standard version of chiral perturbation theory [3] yield numbers in the 
range 

0.2GeV<A3<2GeV. (4) 

The term of order M^ is then very small compared to the one of order M^, so that 
the Gell-Mann-Oakes-Renner formula is obeyed very well. Stern and collaborators 
investigate the more general framework, referred to as "generalized chiral perturbation 
theory", where arbitrarily large values of £3 are considered. The quartic term in eq. (2) 
can then take values comparable to the "leading", quadratic one. If so, the dependence 
of M^ on the quark masses would fail to be approximately linear, even for values of niu 
and niti that are small compared to the intrinsic scale of QCD. A different bookkeeping 
for the terms occurring in the chiral perturbation series is then needed [8] - the standard 
chiral power counting is adequate only if £3 is not too large. 



5. QUARK MASS DEPENDENCE OF M^^ AND F^^ 

The behaviour of the ratio M^/M^ as a function of m is indicated in fig. 1, taken from 
ref. [14]. The fact that the information about the value of A3 is very meagre shows 
up through very large uncertainties. In particular, with A3 ~ 0.5 GeV, the ratio M^/M^ 
would remain close to 1, on the entire interval shown. Note that outside the range (4), 
the dependence of M^ on the quark masses would necessarily exhibit strong curvature. 
The figure illustrates the fact that brute force is not the only way the very small values 
of niu and nij observed in nature can be reached through numerical simulations on a 
lattice. It suffices to equip the strange quark with the physical value of m^ and to measure 
the dependence of the pion mass on mj,, m^ in the region where Mt^ is comparable to Mk- 
A fit to the data based on eq.(2) should provide an extrapolation to the physical quark 
masses that is under good control^ Moreover, the fit would allow a determination of the 
scale A3 on the lattice. This is of considerable interest, because that scale also shows 
up in other contexts, in the KK scattering lengths, for example. For recent work in this 
direction, I refer to [17, 18]. 



^ The logarithmic singularities occurring at next-to-next-to-leading order are also known [15] - for a 
detailed discussion, I refer to [16]. 




FIGURE 1. Dependence of the ratios Fti/F and M^/M^ onm = j{mu + nid). The strange quark mass 
is held fixed at the physical value. The vertical line corresponds to the physical value of m. 



For the pion decay constant, the expansion analogous to eq. (2) reads 



Fn 



i + T^ + oCA'^) 
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M2 



(5) 



In this case, the relevant effective coupling constant is known rather well: chiral symme- 
try implies that it also determines the slope of the scalar form factor of the pion, 

F,(0 = (7r(/)|MM+Jrf|7rb))=F,(0){l + i(r2),? + C>(?2)} . 

As shown in ref. [3], the expansion of (r^) ^ in powers of m„, mj starts with 

6 f. 13 



k^s 



(4;rF)^ 



12 



+ 0(M^ 



(6) 



Analyticity relates the scalar form factor to the 7 = 5-wave phase shift of 7171 scat- 
tering [19]. Evaluating the relevant dispersion relation with the remarkably accurate 
information about the phase shift that follows from the Roy equations [16], one finds 
(r^)^ = 0.61 ± 0.04fm2. Expressed in terms of the scale A4, this amounts to 



A4 = 1.26±0.14GeV. 



(7) 



Fig. 1 shows that this information determines the quark mass dependence of the decay 
constant to within rather narrow limits. The change in F,i occurring if m is increased from 
the physical value to ^ m^ is of the expected size, comparable to the difference between 
Fk and F71. The curvature makes it evident that a linear extrapolation from values of 
order m ~ ^ m^, down to the physical region is meaningless. 



6. 7171 SCATTERING 



The experimental test of the hypothesis that the quark condensate represents the leading 
order parameter relies on the fact that (0| qq |0) not only manifests itself in the depen- 



dence of the pion mass on niu and m^, but also in the low energy properties of the nn 
scattering amplitude. 

At low energies, the scattering amplitude is dominated by the contributions from the 
S- and P-waves, because the angular momentum barrier suppresses the higher partial 
waves. Bose statistics implies that configurations with two pions and i = are symmetric 
in flavour space and thus carry either isospin / = or / = 2, so that there are two distinct 
5-waves. For i = 1, on the other hand, the configuration must be antisymmetric in 
flavour space, so that there is a single P-wave, / = 1 . If the relative momentum tends to 
zero, only the S-waves contribute, through the corresponding scattering lengths Uq and 
Gq (the lower index refers to angular momentum, the upper one to isospin). 

As shown by Roy [20], analyticity, unitarity and crossing symmetry subject the par- 
tial waves to a set of coupled integral equations. These equations involve two subtraction 
constants, which may be identified with the two S-wave scattering lengths «[], ^q. If these 
two constants are given, the Roy equations allow us to calculate the scattering amplitude 
in terms of the imaginary parts above 800 MeV and the available experimental informa- 
tion suffices to evaluate the relevant dispersion integrals, to within small uncertainties 
[21]. In this sense, Uq, Uq represent the essential parameters in low energy nn scattering. 

As a general consequence of the hidden symmetry, Goldstone bosons of zero momen- 
tum cannot interact with one another. Hence the scattering lengths Uq and Uq must vanish 
in the symmetry limit, niu, nid -^ 0. These quantities thus also measure the explicit sym- 
metry breaking generated by the quark masses, like M^. In fact, Weinberg's low energy 
theorem [22] states that, to leading order of the expansion in powers of m„ and nid, the 
scattering lengths are proportional to M^, the factor of proportionality being fixed by the 
pion decay constant:^ 

Chiral symmetry thus provides the missing element: in view of the Roy equations, 
Weinberg's low energy theorem fully determines the low energy behaviour of the Tin 
scattering amplitude. The prediction (8) corresponds to the dot on the left of fig. 2. 

The prediction is of limited accuracy, because it only holds to leading order of the 
expansion in powers of the quark masses. In the meantime, the chiral perturbation series 
of the scattering amplitude has been worked out to two loops [23]. At first nonleading 
order of the expansion in powers of momenta and quark masses, the scattering amplitude 
can be expressed in terms of F^i, M^ and the coupling constants £i, . . . ,£4 that occur in 
the derivative expansion of the effective Lagrangian at order p^. The terms (,\ and £2 
manifest themselves in the energy dependence of the scattering amplitude and can thus 
be determined phenomenologically. As discussed in section 5, the coupling constant £4 
is known rather accurately from the dispersive analysis of the scalar form factor. The 
crucial term is £3 - the range considered for this coupling constant makes the difference 
between standard and generalized chiral perturbation theory. In the standard framework, 



The standard definition of the scattering length corresponds to ao/A^jr- It is not suitable in the present 
context, because it differs from the invariant scattering amplitude at threshold by a kinematic factor that 
diverges in the chiral limit. 
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FIGURE 2. 5-wave scattering lengths. The Roy equations only admit solutions in the "universal band", 
spanned by the two tilted lines. The dot indicates Weinberg's leading order result, while the small ellipse 
includes the higher order corrections, evaluated in the framework of standard chiral perturbation theory. In 
the generalized scenario, there is no prediction for Aq, but there is a correlation between ajj and Aq, shown 
as a narrow strip. The triangle with error bars indicates the phenomenological range permitted by the old 
data, fl|] = 0.26 ± 0.05, ag = -0.028 ±0.012 [27]. 



where the relevant scale is in the range (4), one finds that the leading order result is 
shifted into the small ellipse shown in fig. 2, which corresponds to [24, 25]: 



a'o 



0.220 ±0.005 



al 



-0.0444 ±0.0010 



(9) 



The numerical value quoted includes the higher order corrections (in the standard frame- 
work, the contributions from the corresponding coupling constants are tiny). 

The corrections from the higher order terms in the Gell-Mann-Oakes-Renner relation 
can only be large if the estimate (4) for A3 is totally wrong. As pointed out long ago 
[26], there is a low energy theorem that holds to first nonleading order and relates the 
S'-wave scattering lengths to the scalar radius: 

2'.g-5«-^{l4Mi(.^)., + ^}+0(*3). (.0, 

In this particular combination of scattering lengths, the term £3 drops out. The theorem 
thus correlates the two scattering lengths, independently of the numerical value of A3. 
The correlation holds both in standard and generalized chiral perturbation theory. The 
corrections occurring in eq. (10) at order m^ have also been worked out. These are 
responsible for the fact that the narrow strip, which represents the correlation in fig. 2, is 
slightly curved. 



7. IMPACT OF THE NEW K DECAY DATA 



The final state interaction theorem implies that the phases of the form factors relevant for 
the decay K -^ nnev are determined by those of the 7 = 5-wave and of the P-wave of 



elastic Tin scattering, respectively. Conversely, the analysis of the final state distribution 
observed in this decay yields a measurement of the phase difference d{s) = 5q (s) — 5| (s) , 
in the region 4M^ < s < M\. As discussed above, the Roy equations determine the 
behaviour of the phase shifts in terms of the two 5-wave scattering lengths. Moreover, 
in view of the correlation between the two scattering lengths, a^ is determined by a^, 
so that the phase difference 8(5) can be calculated as a function of ^q and q, where q is 
the cm. momentum in units of Mjt, s = 4M^(1 + q^). In the region of interest (q < \, 
0.18 < «[] < 0.26), the prediction reads 
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(11) 



Z? = 0.2527 + 0. 15 lAa[| + l. 14 (Aa|j)2 + 35.5 (Aa|j 
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c = 0.0063-0. 145 Aa 
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-0.0096 



with Aaj] = «o — 0.22. The uncertainty in this relation mainly stems from the experimen- 
tal input used in the Roy equations and is not sensitive to a^: 



0.0035^7^+0.0015.7^ 
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FIGURE 3. Phase relevant for the decay K -^ nnev. The three bands correspond to the three indicated 
values of the 5-wave scattering length a^. The uncertainties are dominated by those from the experimental 
input used in the Roy equations. The triangles are the data points of ref. [28], while the full circles 
represent the E865 results [13]. 



The prediction (11) is illustrated in fig. 3, where the energy dependence of the phase 
difference is shown for ^q = 0. 1 8, 0.22 and 0.26. The width of the corresponding bands 
indicates the uncertainties, which according to (12) grow in proportion to q^ - in the 
range shown, they amount to less than a third of a degree. 

The figure shows that the data of ref. [28] barely distinguish between the three values 
of Uq shown. The results of the E865 experiment at Brookhaven [13] are significantly 
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FIGURE 4. Ke^ data on the scattering length Oq. The triangles are the data points of ref. [28], while the 
full circles represent the E865 results [13]. The horizontal band indicates the statistical average of the 1 1 
values for Aq shown in the figure. 



more precise, however. The best fit to these data is obtained for «[] = 0.218, with y} = 5.1 
for 5 degrees of freedom. This beautifully confirms the value in eq. (9), obtained on the 
basis of standard chiral perturbation theory. There is a marginal problem only with the 
bin of lowest energy: the corresponding scattering lengths are outside the region where 
the Roy equations admit solutions. In view of the experimental uncertainties attached 
to that point, this discrepancy is without significance: the difference between the central 
experimental value and the prediction amounts to 1.5 standard deviations. Note also that 
the old data are perfectly consistent with the new ones: the overall fit yields a^ = 0.221 
with %^ = 8.3 for 10 degrees of freedom. 

The relation (11) can be inverted, so that each one of the values found for the phase 
difference yields a measurement of the scattering length Aq. The result is shown in fig. 4. 
The experimental errors are remarkably small. It is not unproblematic, however, to treat 
the data collected in the different bins as statistically independent: in the presence of 
correlations, this procedure underestimates the actual uncertainties. Also, since the phase 
difference rapidly rises with the energy, the binning procedure may introduce further 
uncertainties. To account for this, the final result given in ref. [12], 



?[^ = 0.221 ±0.026, 



(13) 



corresponds to the 95% confidence limit - in effect, this amounts to stretching the 
statistical error bar by a factor of two. 

We may translate the result into an estimate for the magnitude of the coupling constant 
£3: the range (13) corresponds to \l^\ < 16. Although this is a coarse estimate, it implies 
that the Gell-lVIann-Oakes-Renner relation does represent a decent approximation: more 
than 94% of the pion mass stems from the first term in the quark mass expansion (2), 
i.e. from the term that originates in the quark condensate. This demonstrates that there 
is no need for a reordering of the chiral perturbation series based on SU(2)« x SU(2)£. In 
that context, the generalized scenario has served its purpose and can now be dismissed. 



A beautiful experiment is under way at CERN [29], which exploits the fact that tt+tt" 
atoms decay into a pair of neutral pions, through the strong transition tt+tt^ ^'kP'kP. 
Since the momentum transfer nearly vanishes, only the scattering lengths are relevant: at 
leading order in isospin breaking, the transition amplitude is proportional to aQ—a^. The 
corrections at next-to-leading order are now also known [30]. Hence a measurement of 
the lifetime of a tt+tt^ atom amounts to a measurement of this combination of scattering 
lengths. At the planned accuracy of 10% for the lifetime, the experiment will yield 
a measurement of the scattering lengths to 5%, thereby subjecting chiral perturbation 
theory to a very sensitive test. 
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